#include <iostream>
#include <stdfin/distribution/gbm_first_hitting_time.hpp>

int main()
{
    for (double t=0, dt=0.01; t<4+dt/2; t+=dt)
        std::cout 
        << t << " " 
        << stdfin::gbm_first_hitting_time_pdf(100.0, 0.05, 0.20, 130.0, t) 
        << std::endl;
    return 0;
}